Passive Retrieval of Rayleigh Waves in Disordered Elastic Media 
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When averaged over sources or disorder, cross-correlation of diffuse fields yield the Green's function 
between two passive sensors. This technique is applied to elastic ultrasonic waves in an open 
scattering slab mimicking seismic waves in the Earth's crust, ft appears that the Rayleigh wave 
reconstruction depends on the scattering properties of the elastic slab. Special attention is paid to 
the specific role of bulk to Rayleigh wave coupling, which may result in unexpected phenomena like 
a persistent time-asymmetry in the diffuse regime. 



I. INTRODUCTION 

Whatever the type of waves involved, knowing 
the Green's function of an heterogeneous medium is 
the key to many essential applications like imaging, 
communication or detection. In the last twenty years 
or so, mesoscopic physics has intensely studied wave 
phenomena in strongly disordered media: weak and 
strong localization, radiative transfer and diffusion 
approximation etc. [lj. But the exact Green's function 
of a complex medium is not easily tractable, and usually 
theoreticians study ensemble-averaged quantities like 
statistical correlations of intensities or wave fields. 
Moreover, from an experimental point of view, it is 
not always possible to measure the Green's functions 
of a complex medium because it requires controllable 
arrays of sources and receivers that do not perturb the 
medium. In the context of laboratory ultrasound, this 
is (nearly) routine; but in other fields of wave physics 
like seismology, distant sources are not controllable. In 
that respect, a good deal of publications followed recent 
works by Weaver. He proposed to cross-correlate the 
diffuse wave fields obtained at two passive sensors and 
showed it yields the elastic Green's function between 
the two receivers, as if one of the receivers was a source 
0, [![■ That correlations performed on passive sensors 
should yield the wave travel times is not that new. This 
principle was applied to helioseismology in the 90 's where 
it provided tomographic images of the Sun's interior 
0, Q. Beyond travel time reconstruction, Weaver's 
experimental retrieval of exact Green's functions was a 
real breakthrough. 

Weaver's work has been followed by different con- 
tributions. Some theoretical works are based on the 
crgodic approximation where time and ensemble average 
coincide. This approach is very useful for estimating the 
role of scattering in the correlation asymmetry [y, 
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When several sources are available, another possibility 
is to average correlations over sources without moving 
the receivers. This was done in underwater acoustics 
p|. Sabra et si @ showed the possibility of recover- 
ing the entire Green's function of the sea waveguide 
(especially the late contributions of multiples) and 
proposed a model for estimating the signal-to-noise 
ratio of the correlations. Derode et al [10] proposed 
to interpret the Green's function reconstruction in 
terms of a Time-Reversal analogy, and showed that 
all the benefits of time-reversal devices in multiple 
scattering media could be fruitfully applied to "passive 
imaging" from correlations [Til . [l3 | . This idea is based 
on the mathematical principle of the representation 
theorem [ljj, equivalently referred to as the Helmholtz- 
Kirchhoff theorem if a perfect series of receivers 

has been sensing the wave field for ages, then one can 
mathematically have access to the wave field anytime 
and anywhere in the area enclosed by the sensors. 
Very recently, this principle was fully developed and 
applied to elastic waves in open media by Wapenaar [l5| . 

From the very start of Weaver's work, the correlation 
principle was successfully applied to seismic waves 
[13, E, Gil- The most energetic part of the Green's 
function between two seismic stations (Rayleigh and 
Love wave trains) was retrieved from passive correlations 
of either coda waves or records of seismic noise. As we 
mentioned earlier, obtaining impulse responses without 
a controllable source is of high interest in seismology 
since it gives the possibility of simulating very energetic 
and punctual earthquakes everywhere around the Earth, 
and do imaging without a controllable source ("passive 
imaging"). Results are quite encouraging although 
several theoretical problems remain unsolved. They 
concern issues very specific to the Earth. In particular 
elastic sources (earthquakes) are naturally never ar- 
ranged in such a way that they would perfectly surround 
a couple of seismic stations. In addition, for distant 
earthquakes, the propagation directions of incident bulk 
waves are mostly vertical, and in a vertically layered 
medium they would not couple with Rayleigh waves. 
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FIG. 1: 10 6 — > 1 scaled analogy between the Earth crust and 
our ultrasonic experiment. In the crust, the scattering mean 
free path was estimated in Mexico [l(| at 1 Hz: I* = 30 km. In 
the aluminum waveguide we numerically found t = 5.5 mm 
at 1 MHz. 



So, why do we observe a Rayleigh wave train in the 
correlation of waves generated by bulk sources? The 
Green's Function retrieval is linked to equipartition, 
which is due to wave scattering. So, what is the influ- 
ence of scattering within the Earth crust in that process? 

Here we present experimental results obtained in the 
lab with laser-induced and laser-detected ultrasonic 
waves propagating in a heterogeneous slab mimicking 
the Earth crust. Lots of previous experimental articles 
applied the "passive imaging" technique to acoustic 
waves. Here we propose to investigate the emergence of 
the Green's function in the correlations of elastic waves 
propagating in a solid heterogeneous medium with a 
free surface. Because field experiments are tedious and 
natural environment arc mostly unknown (especially 
the scattering properties), we propose to build an 
Earth crust model at the scale 1/10 -6 (presented in 
fig. [J). Waves will be sensed at ultrasonic frequencies 
at the free-surface of an elastic open medium. In our 
experiment surface waves are not initially excited. This 
is different from the work of Malcolm et al Q, where 
ultrasonic Rayleigh waves were generated on the same 
surface they were measured. In addition they used 
a finite cylindrical medium with possibly round trip 
wave trains whereas our experiment is conducted in a 
nearly open medium. In our configuration, we would 
not expect the Rayleigh wave to be reconstructed in 
the correlation, except if scattering is present and mode 
conversion between Rayleigh and bulk waves occurs. 
To verify this assumption and study the role of mode 
conversion, we used two samples of identical dimensions, 
one with scatterers the other without. 



The next section describes the experimental setup and 
the propagation medium. Section III presents a short 
theoretical study of the scattering properties (scattering- 
cross sections a of a single scatterer and the transport 
mean- free paths l* P and i* s of the heterogeneous sample) . 
This study is supported by experimental measurements. 
In section IV, the passive imaging technique is applied to 
records acquired at the free surface. Time and frequency 
analysis are proposed, and a brief discussion on the time- 
symmetry of the correlations concludes the article. 



II. EXPERIMENTAL SET-UP 
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FIG. 2: Experimental setup : The propagation medium is 
an aluminum block drilled with 54 vertical cylindrical holes 
(diameter: 4 mm). A Q-switched Nd:YAG laser shoots on 
the bottom side 100 pulses (24 ns duration and 280 MW/cm 2 
intensity each). On the top side a heterodyne interferom- 
eter senses the vertical displacement of the aluminum/air 
interface. This measure is repeated at 7 different locations 
(Xq — mm^ Xq = 60 mm) for each source position. Plas- 
ticine was stuck on the edges and the four lateral sides to 
mimic absorbing boundary conditions and avoid the genera- 
tion of Rayleigh waves by mode conversion at the edges. 

The experimental setup is depicted in fig. [5] It was 
designed to mimic the propagation and scattering of 
elastic waves through the Earth crust, at ultrasonic 
frequencies (0.8-3.2 MHz). A duraluminum slab whose 
dimensions are roughly 10 -6 those of the crust was used. 
Fifty-four cylindrical holes (radius a = 2 mm) were 
drilled at random along direction y so that the waves 
propagating through the slab could undergo multiple 
scattering. The 2-D spatial Fourier transform of the 
hole positions was calculated and is almost perfectly 
flat in the range of ultrasonic wavelengths involved in 
the experiment, which confirms the absence of spatial 
correlation between the holes. The density of scatterers 
was n = 0.0105 mm -2 . Ideally, the slab should have 
had infinite dimensions along x and y. To approach 
this condition, we stuck a thick layer of dense plasticine 
on the lateral sides and edges of the duraluminum 
sample. It was aimed at creating absorbing boundary 



3 



conditions and avoiding the generation of Rayleigh 
waves by mode conversion at the edges. The energy 
decay time r a b s was initially found to be 23,000 /xs. 
With the plasticine it decreased to 120 /is (see next 
section for a detailed discussion on the absorption time). 
We therefore simulated an open slab in the x and y 
direction with a free surface at the top. Rigorously, 
the earth crust is a waveguide that partially leaks 
energy through the Moho to the underlying mantle. 
To perfectly match this feature we should have placed 
a infinite medium with a different impedance at the 
bottom side of the aluminum slab. Yet we think that 
our results and conclusions do not suffer too much 
from this omission: the central point of our set-up 
is to mimic an elastic and scattering medium, infinite 
in the horizontal directions with a free surface on the top. 

The source we employed to simulate earthquakes was 
a Q-switched Nd:YAG laser that shot 24 ns pulses at the 
bottom side of the slab (each pulse energy : 9 mJ). Two 
regimes of elastic wave generation are possible with a 
laser source [2(| • When the surfacic intensity of the pulse 
is weak (< 15 MW/cm 2 ), the surface is locally heated up 
and its dilation creates Rayleigh waves (thermoelastic 
regime). At higher intensity, the laser evaporates part 
of the metal. In this ablation regime, both Rayleigh and 
bulk waves are generated. Theoretical radiation patterns 
are displayed on fig. [3] for compressional and shear waves. 
Rigorously, such a source is not truly reproducible since 
the laser impact can damage the surface. To make the 
experiment as reproducible as possible while staying 
in the ablation regime, the shot intensity was no more 
than 280 MW/cm 2 . A 1 ms record was acquired 100 
times without any observable change. In the following 
experiments, for a satisfactory signal to noise ratio, each 
impulse response was averaged over 100 consecutive 
shots. 

As to the detection of the free surface motion, it 
was achieved with a contactless and quasi-punctual 
device: a heterodyne optical interferometer developed 
by Royer et al [21| which has the advantage of a very 
broadband response (20 kHz-45 MHz) and a sensitivity 
of 10- 4 A/ VlH. It was mounted perpendicularly to the 
slab and then provided us with the absolute vertical 
component of the free surface displacement (top side), 
with a fine spatial resolution; the size of the laser spot 
was ~ 100 /xm whereas the typical elastic wavelengths 
here are ranging between 1 and 10 mm. This is similar to 
seismology, where sensors are nearly punctual compared 
to the wavelengths considered (several kilometers at 
1 Hz). However seismic sensors usually provide time 
records of the three components of the displacement 
field. Here the interferometer only measured the vertical 
movements of the free surface. 

In an elastic body, three different kinds of wave po- 
larization are possible. Compressional (or longitudinal) 



waves are analogous to acoustic waves in fluids (velocity 
vp = 6.32 mm/ixs in duraluminum) . Shear (trans- 
verse) waves have two possible polarizations (velocity 
vs = 3.13 mm//is): one we call SV (Vertical) in the x-z 
plane (see fig. [1]) and one SH (Horizontal) in the x-y 
plane. SH waves have no contribution in the z direction 
and therefore will not be detected by the interferometer. 
In addition to bulk waves, surface waves exist but here 
only Rayleigh waves (vr = 2.9 mm//zs) will be taken 
into account since the others cannot be detected (no 
vertical displacement). The shortest wavelength in the 
aluminium slab is 0.9 mm, which is much greater than 
the duraluminum alloy grain size. Since the orientation 
of the grains is random, we consider the alloy to be 
isotropic for elastic waves in the frequency band of 
interest. Scattering at the grain edge is presumably also 
negligible compared to scattering by the void cylinders. 
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FIG. 3: Directivity pattern (linear scale) of the laser source. 
Rayleigh waves are not taken into account since they are ab- 
sorbed at the edges. The experiments were carried out in 
the ablation regime, where bulk waves are much more ener- 
getic than in the thermoelastic regime. Compressional wave 
directivity: thick line, shear wave directivity: dotted line. 

The overall translational symmetry along y of both 
the free surface and the cylindrical scatterers avoid any 
coupling between SH mode and the other SV and P 
modes. Therefore SH waves will not be considered in 
our article and SV waves will be referred to as S (shear) 
waves. The wave propagation in our experiment will be 
treated as 2-dimensional (and quasi 1-D for the surface 
Rayleigh waves). Waves initially propagating in the y 
direction are rapidly absorbed by the plasticine and lost. 

The laser source and the laser interferometer could 
be translated independently: 35 sources and 7 sensing 
positions were used during the experiment, providing 
us with a set of 35 x 7 impulse responses. A typical 
waveform is depicted on fig. 0J Around 1 MHz it is 
lasting nearly 800 fis, and shows a long diffusive decay 
comparable to the seismic coda. Due to the strong 
scattering on the cylindrical cavities, no top-bottom 
reflection was observed in the data. This confirms the 
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FIG. 4: Top: typical waveform obtained through the scatter- 
ing slab around 1 MHz. The sensitivity of the heterodyne in- 
terferometer is 10 _4 A/V Hz corresponding the optimal level 
of optical reflection on the sensed surface. After averaging 
over 100 records we reached the precision of 10~ 2 A. Bottom: 
intensity averaged over several source/sensor positions. 




FIG. 5: Differential scattering cross-sections of a cylindrical 
cavity calculated for a compressional (P) or shear (S) incident 
plane wave around 1.2 MHz (top) and 2.4 MHz (bottom). The 
scattered wave is either compressional (P) or transverse (S). 
Each pattern is normalized by its maximum. 



high diffusive nature of the propagation in the scattering 
slab. The relevant scattering properties are discussed 
and evaluated in the next section. 



III. WAVE SCATTERING AND TRANSPORT 
PROPERTIES 

A. Scattering cross-section of an empty cylinder 

In order to evaluate the amount of scattering and mode 
conversion, we calculated the differential cross-section 
and the total scattering cross-section a of a cylindrical 
void in a elastic medium excited by a compressional or 
shear plane wave. A brief description of the calculation is 
given in the Appendix. For a detailed derivation we refer 
to [13, HH, [24[ ■ The differential scattering cross-section 
gives the angular distribution of the scattered surfacic 
intensity, normalized by the incident surfacic intensity. 
The total elastic cross-section is a = J ^d9. In 2D it 
has the dimensions of a length. It corresponds to the 
scattering strength of an object at a given frequency. In 
an elastic medium, mode conversion can occur and differ- 
ent cross-sections must be considered. In the case of an 
incident compressional wave, they are noted upp, aps 
and <7p = <7pp + (7 ps, respectively for the P to P, P to 
S and total P elastic cross-sections. We also calculated 
the elastic cross-sections for an incident shear wave (S), 
osp, &SS an d a s = &sp + &ss- 

The differential cross-sections plotted on fig. [5] have 
been computed at 1.2 MHz and 2.4 MHz frequencies. 
The elastic scattering sections are plotted in fig. [5] for 



frequencies ranging from 0.1 MHz to 200 MHz. In 
average in the frequency band of interest (0.8-3.2 MHz), 
we obtained up = 9.2 mm. This value is comparable 
to measurements by White 25]. The same calculations 
were conducted for an incident shear wave, we found an 
average of as — 12 mm in the 0.8-3.2 MHz frequency 
band. 
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FIG. 6: Elastic (er) scattering cross-sections calculated for 
shear (S) and compressional (P) plane wave impinging on 
a cylindrical void with diameter 4 mm. Between 0.8 and 
3.2 MHz, scattering is stronger for shear waves. At high fre- 
quencies the elastic cross-sections tend to the limit of twice 
the geometrical diameter. 
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B. Transport properties 

When the elastic wave propagates through the alu- 
minum slab drilled with holes, it undergoes multiple scat- 
tering. Let ip(t) (resp. fo{t)) be the vertical displace- 
ments sensed at the free surface through the scattering 
(resp. homogeneous) medium. Classically, this field is 
split into two contributions: the coherent and the inco- 
herent part. The coherent wave is the ensemble-averaged 
field (ip{t)) (averaged over disorder configurations, here: 
cylinder's positions). We underline the difference be- 
tween the coherent and the ballistic wave (i.e. the first 
arrival). For a detailed discussion about scattering ef- 
fect on coherent and ballistic waves, see [2g|. Away from 
resonances, the coherent wave can be roughly thought of 
as an attenuated version of the direct wavefront (po(t). 
When there is no intrinsic dissipation, the energy of the 
coherent wave decays with the slab thickness H as e~ H / e , 
where £ is the elastic mean free path. Assuming a dilute 
set of scatterers, the elastic mean- free path is simply re- 
lated to the scatterers density n and their elastic cross- 
section a: 




From the theoretical scattering cross section calculated 
above, we find £p = 10.3 mm. In order to measure the 
mean-free path experimentally, we used two aluminum 
slabs of exactly the same dimensions. The first served 
as a reference and provided measurements of tpo{t) for 
different source-sensors positions. The second one was 
drilled with holes. By translating the source-receiver de- 
vice along the slab, we achieved something very similar 



to a configurational averaging and measured the energy 
of the coherent wave (ip(t))' 2 . Between 0.8 and 3.2 MHz, 
we obtained £p = 9 ± 2.5 mm from these experiments. 

The intensity of the incoherent part was also studied. 
The time evolution of the averaged incoherent intensity 
I(t) = ((,c(t) 2 ) is governed by another parameter: the 
transport mean free path £* . In an elastic body, transport 
quantities have been theoretically defined by [27[: 

p* = i a s - ggs + cr* PS m 

n {up - a PP ){as - a ss ) - cr PS a SP 



I* = I <?p-<?pp + <?Ip (2) 

n [a P -<Jpp){(Js - <J SS ) -crps a SP 

with a* = J ^§cos{6)d9. It was evaluated numerically: 
£* P ~ £* s = 5.5 mm. In an experiment, this parameter 
is very hard to measure with a reasonable prec ision. 
The coherent backscattering effect [H H, HJ 1H| (also 
referred to as weak localization) does give a direct 
estimation of the transport mean free path but our 
experimental configuration did not allow this special 
measurement since we could not place a laser sensor in 
the vicinity of the laser source. Yet we checked that the 
experimental intensity decay I{t) gives an order of mag- 
nitude for £* that is consistent with the theoretical value. 

For the sake of simplicity we propose a 2-D scalar wave 
model for I(t) [32| . under the diffusion approximation. In 
an infinite slab of thickness H with perfect reflections on 
both sides, the averaged transmitted intensity reads: 



I(X,t) = I < == + Y V ' e «» ™. 

I 



with T a bs the absorption time (taking into account the in- 
trinsic absorption in the aluminum and the lateral leaking 
due to the plasticine) and D the diffusion constant. X 
is the lateral distance between source and receiver. This 
formula is obtained using a modal decomposition of the 
diffusion equation in the z direction. The intensity I(t) 
is a mix of compressional and shear waves, each mode 
traveling with its own parameters (velocity, £* , diffusion 
constant) and interchanging their energy through scat- 
tering events. In our experiment £* s and l* P are of the 
same order. The diffusion constant was approximated 
b y D = 2(1+4/4) {Vptp + 2vs ^ H s) « 40 mm 2 / M s. 
This assumption is valid after a couple of mean free times, 
when the equipartition regime [33j is set. Equiparti- 
tion means that the density of compressional and trans- 



verse modes equilibrates. Considering the specific veloc- 
ities of each mode [13], we infer that 80% of the energy 
is transported by S waves, and only 19% by P waves 
(and an additional 1% for surface waves). Hence our 
best fit (fig. [7]) of the intensity decay in the coda gives 
r abs = 120 ^s±10% and I* = 5 - 20 mm. 



IV. TWO-POINT CORRELATION OF DIFFUSE 
FIELDS 

In this section we focus on the experimental recon- 
struction of the direct Green's function from "passive" 
correlations. The main idea is to correlate diffuse fields 
sensed at two different locations on the top side when 
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FIG. 7: Averaged intensity I(t) and theoretical fit for scalar 
wave diffusion in a 2-D semi-infinite slab with thickness H = 
24 mm. i* — 5.5 mm and r a ba ~ 120 /is. 



a source generates bulk waves at the bottom. Since 
we record the vertical component of the surface dis- 
placements, the two-point correlations should simulate 
a vertical source at the surface, which mainly generates 
surface waves. Indeed, the experimental correlations 
we obtained reveal a wave packet that travels at the 
speed of a Rayleigh wave. We insist that our sources do 
not generate surface waves at the top side of the slab. 
Moreover if a surface wave happened to be generated 
anywhere, it would be completely absorbed by the 
plasticine. Under these conditions no Rayleigh wave 
should travel on the top surface, and no Rayleigh wave 
should be passively retrieved by correlations. Why then 
should passive imaging give rise to a Rayleigh wave train 
in our experiment? 

We propose first to examine the role of scatterers for 
the emergence of the direct Rayleigh wavefront in the 
correlations. To that end we separately correlated coda 
records obtained through two different aluminum slabs: 
the first drilled with holes, the second without. Each im- 
pulse response was lasting w 800 (is before reaching the 
noise level (see fig. H}. We underline that these record 
lengths are far from the Heisenberg time (break time) 
at which the modes of the aluminum block would be re- 
solved (here Tjj w 10 6 /zs) and correlations would natu- 
rally converge to the Green's function. This modal ap- 
proach is unrelevant to our experiment. The records were 
correlated and averaged over the 35 available sources. For 
the scattering slab, this reads: 

35 „t=60Q fis 

(Cij(r)) = / <p{S,Xi,t)<p(S,X jt t + T)dt 

J t—0 [is 

where Xj and Xj are the sensors points (running from 
Xq = mm to Xq = 60 mm along the array). And for 
the homogeneous slab: 



35 ,,(=600 fia 

( C U T )) = E / cpo(S,Xi,t)cpo(S,Xj,t + T)dt 

J t— }±s 
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FIG. 8: Green's function reconstruction for different pairs of 
receivers. Correlations are averaged over the 35 sources and 
filtered in the 0.8-1.6 MHz frequency band. On the left (a) 
are displayed the 7 cross-correlations {Cij(r)} in the diffusive 
aluminum plate. On the right (b) the 7 cross-correlations 
(Cfj(r)} are obtained in an equivalent aluminum block with- 
out any hole, where the mode conversion possibly occurring 
at the edges was suppressed by the plasticine. The insets 
show the summation of the 6 wavefronts from Xo — 10 mm 
to Xq = 60 mm after a time-reduction based on the Rayleigh 
wave speed (dotted line, vr = 2.9 mm//js). 
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FIG. 9: Same figure as fig. [8] except correlations are filtered 
in the 1.6-3.2 MHz frequency band. The signal-to-noise ratio 
in (a) is increased compared to the results obtained at lower 
frequencies. 

To enhance the signal-to-noise ratio, each correlation 
is time symmetrized (— (C(+r)) + (C(— r)}) and 



normalized by its maximum. Results are displayed in 
fig. [5] and [HI A propagating wavefront (traveling at 
the Rayleigh wave velocity Vr = 2.9 mm//is) is clearly 
visible in the presence of scatterers, whereas it does not 
appear in the homogeneous slab. We also summed the 
6 normalized propagating peaks after having delayed 
each signal according to the Rayleigh wave travel time. 
The summation is displayed in the enclosed box on 
each figure. In the scattering slab, its amplitude nearly 
corresponds to the coherent addition of 6 pulses. In the 
homogeneous device, the amplitude of the summation 
is sa 2.5 (incoherent addition of 6 uncorrelated fields). 
We conclude on the necessity of mode coupling due 
to scattering for the Rayleigh wave to emerge from 
the passive correlations of diffuse fields generated by 
bulk waves sources. This is especially relevant for 
applications to seismology. Therefore, it appears once 
again [ll[ that the role of scattering is crucial in " passive 
imaging" . Firstly, because of multiple scattering, at late 
times the equipartition regime can be attained whatever 
the sources/receivers positions. Secondly, because of 
mode conversions due to the scatterers, a Rayleigh wave 
emerges from the passive correlations even though no 
Rayleigh wave was generated by the sources. Note that 
the bandwidth in the upper band record is a little wider 
than in the lower band. This was done to compensate 
for the coda shortening (< 600 [is) so that the product 
TAf is kept constant. 

We can go a little further and catch a glimpse of 
the time symmetry properties. We performed the 
correlations into two consecutive time windows: from 
to 45 [is and from 45 [is to 600 [is, and did not 
time-symmetrize the correlations (fig. 10 and 11). 45 [is 
is twice the time after which the diffuse energy spreads 
homogeneously along the array of receivers (length 
L = 30 mm) : L 2 /AD w 22 [is in an open 2D scattering 
medium. The time series were filtered in two frequency 
bands: 0.8-1.6 MHz and 1.6-3.2 MHz. 

In the first time-window, from to 45 /us, correlations 
are asymmetric in time. This means that the causal part 
(t > 0) and the acausal part (r < 0) of the correlations 
are different (see left part of fig. [TU] and fig. [TTj) . In 
the causal part, a Rayleigh wavefront is clearly visible 
whereas noise is dominating the acausal part. This 
is due to the preferential direction of Rayleigh wave 
propagation (waves traveling from Xq to A3 in our 
experiment). There is a net flux of energy from Xq to 
Xi, X 2 and X 3 (distances 10, 20 and 30 mm in fig. [TU1 
and fig. [TTj) . This flux is due to the uneven distribu- 
tion of sources in comparison to the receiver couples 
X — Xi, X — X 2 and X — A 3 : most of the sources 
are on the Xo's side. At the early times of the coda 
(from to 45 [is) , the diffusion regime is not yet attained. 

Later in the coda, from 45 [is to 600 [is, the wave 
field in the bulk of the aluminum slab is very likely 
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FIG. 10: Evolution of the asymmetry in the reconstructed 
Green's function in the low frequency regime (0.8-1.6 MHz), 
for early (a) and late (b) times. At early times, the waves 
are exciting scatterers on one side of the receiving array more 
than the other. The energy flux is clearly going from the 
source to the rest of the medium. After 45 [is the field is 
equipartitioned and the scattering halo fills the array, leading 
to a more symmetric correlation. 



to be equipartitioned. In the low frequency band 
(0.8-1.6 MHz), the time-symmetry of the correlation 
is indeed restored [H, [36| (see right part of fig. 1X0]) : 
Rayleigh waves travel in all directions. Nevertheless 
and surprisingly, the asymmetry persists in the high 
frequency band (from 1.6 MHz to 3.2 MHz, see right 
part of fig. Ill[) . To interpret this observation, we have 
carefully studied the location of the scattering sources 
around the array. In the high frequency band, the 
Rayleigh wavelength is ~ 1.5 mm. The generation of 
Rayleigh waves by scattering necessarily occurs in the 
first half- wavelength beneath the free surface [35[ . In 
our scattering slab, one hole was nearly showing on the 
surface (position X < 0), another one was 0.61 mm 
beneath (position X > 60 mm) , the others being located 
much deeper. Rayleigh wave trains are mainly generated 
by the hole nearest to the surface, then propagate along 
the array of receivers (from Xq to Xq point). These 
waves are almost not perturbed (attenuated) until they 
reach the edges of the slab and the absorbing plasticine. 
They contribute to a very clear propagating pulse in 
the positive part of the correlation. The weak coupling 
due to the deeper hole (z = —0.61 mm) on the Xq side 
contributes to a smaller pulse propagating from Xq to 
Xq in the negative part of the correlations. 

This interpretation is in agreement with observations 
in the low frequency band (0.8-1.6 MHz), where the 
average Rayleigh wavelength is 3 mm (fig. [TO]) . At least 
5 holes are present in the first half wavelength and 
should cause significant scattering of Rayleigh waves 
and coupling between surface and bulk waves. This 
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time, the holes are evenly distributed along the sensor 
array. In the late coda, correlations X x Xi, X a x X 2 
and Xq x X3 are nearly symmetric. The time sym- 
metry is obtained thanks to scattering by a symmetric 
distribution of scatterers. Under these conditions, a 
global equipartition among bulk and surface waves is 
guaranteed. Furthermore the reconstructed surface wave 
is strongly attenuated along its path because it senses 
many scatterrers along the array. In addition, these 
scatterers contribute signals around t = in the corre- 
lations, which degrade the reconstruction. We think this 
interpretation explains why the symmetric wavefront is 
much more noisy at low frequencies (fig. [8]), where lots 
of cavities are encountered in the first half wavelength, 
than at higher frequencies where the Rayleigh wavefront 
can propagate freely (fig. fTU|) . 

Finally, we comment on the possible misidentification 
of the waves that are reconstructed in our experi- 
ment. Indeed, if many scatterers are present within a 
wavelength, the overall wave velocity may be different 
(effective medium). In the heterogeneous plate, a 
reduced-speed shear wave might propagate with the 
same wavespeed as a Rayleigh wave in the bulk alu- 
minium plate (i.e. without cavities). In our experiment, 
the hole interspacing is 10 mm on average, which is 
larger than the largest ultrasonic wavelength. Thus, it 
is reasonable to assume that the shear waves propagate 
in the aluminium plate with the same wavespeed as in 
the bulk. The measured wave velocity is 2.9 mm//zs at 
all frequencies (from 0.8 MHz to 3.2 MHz) and indeed 
corresponds to the velocity of the Rayleigh wave. 
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FIG. 11: Evolution of the asymmetry in the reconstructed 
Green's function in the high frequency regime (1.6-3.2 MHz), 
for early (a) and late (b) time windows. At early times, the 
waves are exciting scatterers preferably on one side of the 
receiving array. After 45 /j,s, bulk waves are expected to be 
equipartitioned, but not Rayleigh waves. Mode conversion to 
Rayleigh waves mainly occurs at one scattering cavity, leading 
to an anisotropic energy flux, and an asymmetric correlation. 



V. CONCLUSION 

In this paper were presented laboratory experiments 
of elastic wave propagation in heterogeneous media at 
ultrasonic frequencies. An aluminum slab was made 
quasi-infinite by the use of absorbing boundaries to 
mimic the Earth's crust, in which scattering was ob- 
tained drilling cylindrical cavities. A relatively simple 
theoretical model for wave scattering properties was 
proposed. Wave field generation and detection was 
achieved using contactless and quasi-punctual laser 
devices. The cross-correlation of diffuse fields was 
performed, allowing us to retrieve passively the Rayleigh 
wave between two sensors only when scattering was 
present. Without scatterers, and in the case of bulk 
wave generation and surface detection, no Rayleigh wave 
was reconstructed. This illustrates the role of scattering 
and mode conversion in the Green's function passive 
reconstruction. 

Analysis for different time-windows and frequency 
bands were conducted. Previous works in acoustics [37| 
and seismology 17] observed that the Rayleigh wave re- 
construction was harder with increasing frequency. Here 
we found that the Rayleigh wave reconstruction was more 
efficient in the high frequency band. In addition we ob- 
served that even in the late coda where waves are ex- 
pected to be equipartitioned, asymmetry in the correla- 
tions may remain. Both observations are due to the very 
specific coupling between bulk and Rayleigh waves, which 
occurs if scatterers are present in the first wavelength 
beneath the free surface. We emphasize that equiparti- 
tion of bulk waves does not always mean equipartition of 
surface waves. On the one hand, our experiments show 
the need of scattering to passively retrieve the impulse 
response between two sensors, on the other hand they 
show that scattering occurring between the sensors de- 
grade the reconstructed Rayleigh waves. The trade-off 
between the two effects should be further investigated. 
Though our experiment was designed for seismological 
applications, results should be applicable to other fields 
of wave physics where both surface and bulk waves are 
present. 



APPENDIX A: CALCULATION OF THE 
SCATTERING CROSS-SECTION OF A 
CYLINDRICAL CAVITY. 



Here follows a brief calculation of the wave field scat- 
tered by a cylindrical cavity insonified by a plane com- 
pressional wave. Three displacement potentials are rele- 
vant: <F; is the displacement potential of the incident P 
wave, $ s is the scattered P wave potential and ^f s is the 
scattered S wave potential. They can be expanded as: 
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up = ap^p + ap^s 



$;(r,i) = e m i m J m (k P r)cos{ 



rat) )e 



The corresponding differential scattering cross-sections 
are given by: 



E 

m=0 



A m JT^ 1 )(A;pr)c OS (TO0) e - i " t 



^ m— 



¥.(r,t) - ^ B m iJ^)(fc s r) S m(m^)e- Mt 

where J m and 7?™'' are respectively the Bessel and 
Hankel functions both of first kind and of order m and 
e m is the Neumann factor. Taking into account the null 
traction condition at the surface of the cylinder allows 
the calculation of the A m and B m coefficients. Those co- 
efficients are given in [22|, [2^, [24j . The scattering cross- 
sections are 



2 

kp 



<j P ^P = ^l2\A Q \ 2 +J2\A m \ 2 ] 



= £-1 £ B m i- m sin(me)\ 2 

^ m— 1 
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